Curvature, hybridization, and STM images of carbon nanotubes 
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The curvature effects in carbon nanotubes are studied analytically as a function of chirality. The 
7r-orbitals are found to be significantly rehybridized in all tubes, so that they are never normal to 
the tubes' surface. This results in a curvature induced gap in the electronic band-structure, which 
turns out to be larger than previous estimates. The tilting of the 7r-orbitals should be observable 
by atomic resolution scanning tunneling microscopy measurements. 

PACS numbers: 61.48. +c, 61.16.Ch, 71.20.Tx, 73.61.Wp 



The electronic band-structure of carbon nanotubes has 
been a topic of intense investigation ever since their dis- 
covery in 1991 fl. The basic electronic properties were 
quickly understood by numerical studies of the graphite 
tight binding band structure together with a simple zone 
folding model In graphite the four outer electrons 

of carbon form three sp 2 -hybridized a-bonds and one ir- 
orbital, which gives the conduction band with six Fermi 
points and a linear dispersion around each of them jjj. 
The electronic structure of the nanotube is then deter- 
mined by the chiral wrapping vector along the direction 
(n, m) since this determines whether the Fermi points 
satisfy the nanotube's circumferential boundary condi- 
tions. In that model, tubes with a chiral vector that 
satisfies mod[(n — m)/3] — have their Fermi point in 
the allowed fc-space and thus are considered to be metal- 
lic, while all other tubes are semiconducting pl[3]|. But 
even the "metallic" tubes may open a small gap if the 
bond symmetry is broken due to curvature [^|(|, which 
has been analyzed in analytical studies in terms of a one- 
orbital tight binding approximation 

One difficulty in predicting the effect of curvature on 
the electronic properties has been to determine the ex- 
act bond energies in the curved graphite sheet to arrive 
at an analytical formula. So far it has explicitly been 
assumed that the 7r-orbitals are orthogonal to the tubes 
surface [Q, which is a common overly simplified picture 
that has also been used in a previous report by the au- 
thors We show in this Letter, however, that the tt- 
orbitals are never orthogonal to the surface, and instead 
are rehybridized due to the effect of the lower lying a- 
bonds. Typically such a mixing effect is always expected, 
but early studies have estimated that this band-mixing 
can be neglected ||. However, for the case of metallic 
tubes we find that this mixing plays an important role, 
which is crucial in determining an analytic formula for 
the curvature induced bandgap as a function of chirality 
and curvature. Moreover, we can predict the explicit an- 
gles of the 7r-orbitals relative to the tubes surface, which 
should be observable in atomic resolution pictures from 



Scanning Tunneling Microscopy (STM) experiments. 

Our starting point is the well-defined geometrical 
structure of the carbon nanotubes by describing it in 
terms of an "unrolled" graphite sheet. The <T-bonds lie 
along three vectors which can be expressed in a coordi- 
nate system of the circumferential and translational axes 
(c, t) in terms of the chiral indices (n, m). 
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2. 49 A is the length of the honeycomb unit 
,m + m 2 is the circumference in 



where a 
vector and Ch 
units of a. 

In the regular sp 2 hybridization of the unrolled 
graphite sheet the four atomic wave functions can be 
written as 
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where \s) stands for the atomic s-orbital and \t), |c), \ z) 
denote the p-orbitals along the translational, circumfer- 
ential and normal directions in a nanotube, respectively. 
Here /3j are the angles of the bonds relative to the circum- 
ferential direction (cos /3j = c-Ri). Each carbon atom has 
its own local coordinate system, where the z-direction is 
given by the normal direction to the graphite surface. In 
a nanotube neighboring atoms have a relative angle 2oti 
between their z-directions as shown in Fig. ^ with 
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where r = ac^ /2tt is the radius of the nanotube. We 
call this a geometrical tilting of the 7r-orbitals, which is 
known to induce a curvature gap |?],|8| and is predicted 
to cause a stretching around the circumference of STM 
images M. 
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FIG. 1. Schematic cross-section of a bond Ri in a nan- 
otube with radius r. The 7r-orbitals are no longer normal to 
the tubes surface, but are titled by the hybridization angle 
S = a/4:\/3r (small arrows). 



However, in addition we find an equally important con- 
tribution to the tilting from hybridization, which we will 
discuss next. This hybridization comes from the fact 
that in a nanotube the three cr-bonds are not in the 
same plane, but instead directed towards the positions 
of the nearest neighboring carbon atoms, ie. they are 
tilted down by the angles a, relative to the tangential 
c-direction as shown in Fig. |l|. The hybridization of the 
(T-bonds is therefore changed from the uncurved expres- 
sion for |of ) in Eq. (0) to 

ki> = S t \8) (4) 

+ yl — sf (smPi\t) + cosaicos/3i|c) — sina^ cos/3i|z) ) , 

where the mixing parameters s, (expanded around Sj ~ 
•v/l/3) can be determined by the three orthonormality 
conditions between the a bonds {ci\(Jj) = 5%j. The hy- 
bridized 7r-orbital can now be calculated in terms of the 
local basis of atomic orbitals by using the orthonormality 
conditions 

Moi) = 0. (5) 

In what follows we will only work to lowest order in the 
curvature parameter a/r = 2ir/ch because energy band 
repulsion will give higher order corrections for narrow 
nanotubes. In particular, it has been shown by ab initio 
LDA calculations that the c-bands are pushed up above 
the Fermi points for very small radii r < 2. 4 A [fLOf. The 




FIG. 2. The three cr-bonds around a carbon atom in the 
rolled up graphite lattice. The unhybridized 7r-orbital is 
marked by the black dot in the center as seen from above. 
The predicted hybridization in Eq. (^|) of the 7r-orbital is in- 
dicated by the triple arrow with 5 = a/4\/3r. The example 
shown here corresponds to a (9,3) nanotube seen from a height 
of h « 3a. Note, that the angle 9 in Eq. ^ can in fact be 
chosen to be any of the three angles 0i of the bond vectors 
relative to the c-direction to get the correct hybridization. 

straightforward hybridization analysis here is therefore 
only correct as a lowest order approximation in a/r <C 
1, but gives a useful physical picture for most observed 
nanotubes. 

From Eqs. @ and (§) we can now find the correct 
expression for the 7r-orbitals to lowest order in a/r 

Itt) ~ \z) + — %=- (V2\s) ± sin 30|t) ±cos36'|c) N ) , (6) 
4V3r v / 

where 9 is the so-called chiral 
angle = min(|/3i|, |/3 2 |, |/3 3 |), < 9 < tt/6 as shown 
in Fig. H The different ± signs in Eq. (||) refer to neigh- 
boring atoms A and B in the bipartite graphite lattice, 
since their bonds and the corresponding angles (3i point 
in opposite directions as indicated in Fig. ^. 

The physical interpretation of Eq. (j^) is very intuitive. 
The 7r-orbital is always inclined by the hybridization an- 
gle S of size 5 = a/AySr relative to the normal direction 
to the tube's surface. However, the direction of this in- 
clination rotates with 39 relative to the c-direction as in- 
dicated in Fig. |^. It can be verified that the three angles 
between the 7r-orbital and each of the cr-bonds in Eq. (Q) 
are equal as should be expected by symmetry. Some- 
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FIG. 3. Two neighboring carbon atoms with with the pre- 
dicted hybridization angle S = a/4y/3r and geometrical curva- 
ture effects sin on = 28 cos /% , seen from above. The example 
shown corresponds to the bond Ri in a (9,3) nanotube. 



times it is useful to express the hybridization in Eq. (g) 
in terms of the chiral indices (n, m) using 

sin36> = (n - m)(2m 2 + 5mn + 2n 2 )/2c\ (7) 
cos 39 — 3V3mn(m + n)/2c\ 

At this point we can proceed to calculate the hopping 
matrix elements 7j between neighboring 7r-orbitals which 
will determine the one-orbital tight binding band struc- 
ture. Following the Slater-Koster scheme |12| for calcu- 



lating the matrix elements 7i between tilted 7r-orbitals we 
need to use a common coordinate basis for the two neigh- 
boring atoms. Both the hybridization effect in Eq. (^) and 
the geometrical tilting with ai in Eq. (||) can easily be ex- 
pressed in the coordinate system where the bond vector 
Ri between two neighboring atoms A and B defines the 
x-direction as shown in Fig. || 



\ir) ~ \z) + y25\s) ± (<5cos2/3i — sin ca cos (3i 
± (<5sin2/3i + sin ai sin/%) \y), 



(8) 



where S — a/4^/3r and the angles /% are now defined in 
respect to the A site. Using the notation in Ref. ||, we 
can use the overlap integrals between neighboring s- and 
p-orbitals V ssa ,V spa ,V ppa , VJ, p7r to calculate the hoping 
matrix elements "fi 



^ = V t 
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(9) 



([3 + 8 sin 2 2(3i}V ppn - 2V ssa - 2V2V spa + V ppa 



Here we have also used the second order term for the \z) 
orbital, which contributes to this expression. 

The most interesting aspect of the electronic structure 
in metallic tubes is the size of the gap due to curvature 
as a function of the chiral wrapping vector, which we can 
now calculate directly. After having determined the re- 
hybridized orbitals and hoping integrals, we can typically 
ignore the lower lying bands in further calculations and 
use a simple one-orbital tight binding approximation. In 
this model the gap can be calculated from the positions of 



the Fermi points fci? in the curved graphite sheet which in 
turn are determined in terms of the new hopping matrix 
elements 7^ 0] 
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Using the linear dispersion relation and the position of 
the quantization lines it is then straightforward to derive 
the gap equation [||,|ll| assuming that the tube is metallic 
mod[(n — m)/3] = 



E„ 



2^3 



a 2 
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This surprisingly simple formula reconfirms again the no- 
tion that armchair tubes (n = m, 8 — 0) do not have a 
gap from curvature, while zigzag-tubes (m = 0, 9 = ir/6) 
and all other metallic tubes acquire a gap of order 1/r 2 
[||. For semiconducting tubes with an intrinsic gap of 
order 1/r this geometrical and hybridization correction 
to the gap can typically be neglected, but for metallic 
tubes it is essential to take the proper hybridization into 
account. While the dependence of the gap on the chi- 
ral angle agrees with previous analytic studies [0,§], the 
expression in Eq. (|l^) is by a factor of four larger than 
those estimates, which have not considered hybridization 
effects. 

We now discuss how the rehybridization may be ob- 
servable in STM experiments. Already without rehy- 
bridization the directions of the 7r-orbitals in a curved 
geometry can affect the STM images as predicted in 
Ref. || . In that study only the geometrical effects were 
taken into account, which resulted in an effective stretch- 
ing of the STM image along the circumference by the 
amount x c — > x c (l + h/r), where h is the height of the 
STM tip relative to the nanotube surface. However, the 
tilting S caused by hybridization is of similar size and will 
results in an additional distortion of the STM image in 
both transverse and circumferential directions, but alter- 
nating for the inequivalent A and B atoms 

/ h\ ha 
x c — ► x c 1 H — ± — T=~ cos 60 

x t — > x t ± ^° sin36>, (12) 
4V3r 

where x c and xt measure distances in the circumferential 
and transverse directions, respectively. Interestingly, the 
distortion from the hybridization depends again to linear 
order on h/r in all tubes, as opposed to the curvature gap 
which is a correction of second order and only matters for 
metallic tubes. If we consider this distortion for the case 
of a zigzag tube with 9 = 7r/6, we see that the hexagons 
of the graphite sheet will appear compressed in the trans- 
verse direction, but stretched around the circumference 
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but the results coincide with analytic studies |t||| that 
did not take any hybridization into account. Higher or- 
der effects and energy band repulsion will modify the size 
of the hybridization angle S and therefore also the energy 
gap. However, the direction of the hybridization given by 
36 relative to the circumference is correct as can be seen 
by symmetry arguments. The exact size of the hybridiza- 
tion angle can most reliably be found by analyzing STM 
images as outlined above, which in turn would lead to a 
more reliable estimate of the curvature gap. 




Zigzag Armchair 



FIG. 4. The predicted distortion of an STM image (dashed 
lines). For a zigzag tube we find a pronounced squaring of 
the hexagons and a change of the lattice constants along the 
transverse directions. For the armchair tube the image is only 
affected along the circumferential direction. 

as shown in Fig. || For the armchair tubes 9 = on the 
other hand there is no deformation along the transverse 
direction, but the hexagons are still deformed along the 
circumferential direction as shown in Fig. ^. To observe 
the hybridization effects it would therefore be most ad- 
vantageous to scan along the ridge of a zigzag tube and 
average or Fourier transform the image over a distance 
of several hundred carbon sites, preferably for more than 
one value of the height h. 

In summary we have shown that it can be expected 
that the curvature of carbon nanotubes will result in a 
significant rehybridization of the 7r-orbitals. This rehy- 
bridization will affect the energy gap and will also man- 
ifest itself in a well-defined distortion of atomic resolu- 
tion STM images given by Eq. (||). To lowest order 
in the curvature parameter a/r the hybridization angle 
was found to be 6 = a/Ay/Sr, resulting in an energy 
gap which is significintly higher than previous analyti- 
cal studies |t]|| as well as numerical estimates ||. It 
is not clear to us what assumptions were used to model 
the curvature in the tight binding calculations of Ref. || , 
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